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Abstract 

The contrasting melting behavior of different surface orientations in metals 

can be explained in terms of a repulsive or attractive effective interaction 

between the solid-liquid and the liquid-vapor interface. We show how a crucial 

part of this interaction originates from the layering effects near the liquid 

metal surface. Its sign depends on the relative tuning of layering oscillations to 

the crystal interplanar spacing, thus explaining the orientational dependence. 

Molecular dynamics recrystallization simulations of Au surfaces provide direct 

and quantitative evidence of this phenomenon. 
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Surface melting (SM), i.e. complete wetting of a crystal surface by a thin film of its own 
melt, is very common in nature. For example, all faces of a real system, like rare gas solid 
Ar as well as of its popular model, the fee Lennard- Jones (LJ) crystal 0, are believed to 
melt completely. On the other hand, most metals have at least one close-packed face which 
does not melt, and is crystalline right up to [Q. This non-melting (NM) behavior is often 
explained phenomenologically, in terms of a low solid-vapor (SV) surface free energy 7sv, 
which in this case can be lower than the sum of solid-liquid (SL) plus liquid-vapor (LV) 
interface free energies 7sl + 7lv- Open faces, by contrast, tend to have a higher energy 
when solid and this favors SM, in agreement with observations. At a more microscopic level, 
Chernov and Mikheev made the interesting additional remark that surface non-melting 
could in fact be ascribed to large minima in a strong oscillatory effective interaction between 
the SL and the LV interface. This oscillatory interaction is caused by a layering effect: the 
solid (similarly to a hard wall) induces in the liquid film a layer-like density oscillation 
with periodicity lixjQo [Qo is the wavevector where the liquid structure factor S{ff) has its 
strongest peak], which in turn generates the interaction. Binding of the two interfaces in 
the deepest oscillation at £ = is then responsible for a minimum of 7sv and leads to NM. 
Alternatively, if the oscillations are weak, they can be washed out by capillary waves, the 
two interfaces unbind and there will be SM. 

The purpose of this Letter is to argue theoretically, and to demonstrate through detailed 
simulations, that in non-LJ systems, particularly in metals, there are in reality not one, but 
two distinct layering oscillations — one with periodicity 271 /Qo, tied to the LV interface, and 
one with the periodicity of the interlayer spacing a, tied to the SL interface — which overlap 
and interfere inside a liquid metal film. Tuning and detuning of the two oscillations depends 
on orientation, and may dramatically change the strength of the interaction, thus affecting 
deeply the surface melting behavior. In particular, while as expected well-tuned layering 
effects disfavor surface melting (and can make the close-packed surface very resistant to 
overheating), the same does not happen for large detuning, where SM is favored. 

We will proceed in three stages. First, we give a qualitative, analytical argument show- 
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ing that layering effects are strongly dependent on the orientation. Second, we make this 
more concrete through a model ID wetting calculation adapted to our case. Lastly, we 
present detailed surface recrystallization molecular dynamics (MD) non-equilibrium simula- 
tions which illustrate realistically the effect of attractive layering forces for a fcc(lll) NM 
face, in contrast with a fcc(llO) SM face where repulsion dominates. 

Consider a liquid film of thickness £, wetting its own solid. The effective SL-LV in- 
terface interaction free energy can be written as f{£) = /sr(^) + /l(^) where /sr is a 
short-range effective interaction, accounting for the merging of the two interfaces at close 
contact, and /l is a residual part due to layering effects, which can be characterized as fol- 
lows. The underlying solid propagates into the liquid film a first density oscillation roughly 
of the form Spsh{z) ~ ksexp[—{z — zsL)/^s]cos{2nz/a) which has the periodicity of the 
(face-dependent) interlayer spacing a, and decays outwards from the SL interface center 
zsL into the liquid {z grows going from solid to liquid to vapor) within some characteristic 
length 0- The LV interface of a metal originates however a second density oscillation 
6puv{z) ^ kyexp [{z — 2;lv)/^v] cos [{z — zlv) Qo] which propagates inwards, with the typi- 
cal liquid periodicity, 2-71 /Qo, and with a decay length related to the peak width 6Q in 
the liquid structure factor S{Q), 2tt/^y ~ 6Q 0. This second oscillation is not universally 
present: it is only expected for a strongly relaxed and contracted surface, like that of a 
metal, which represents a heavy disturbance for the liquid below 0. In this sense, the LV 
interface in a metal is similar to a hard wall. By contrast, at the more disordered LJ liquid 
surface this oscillation (still present in principle) is below 2% and practically irrelevant 
at Tm- A quantitative picture obtained by simulation (fig. ||, described later) shows that a 
large oscillation is present for Au and absent in LJ. When £ = zu/ — zsl decreases, the two 
density oscillations overlap, giving rise, within linear response, to 

fdi) = \j d^rd'r'Civ - r')5psL(r)5pLv(r') (1) 

where the integration should be restricted to the liquid region z^l < z < z^sj, and C(r — r') 
is the liquid density response function. Taking as a simple approximation for large r the 
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form C{r) = Co exp(— r/^) cos(Qot), the integral can be carried out explicitly. The (111) 
surface corresponds to Qoa = 27r, the (110) surface to Qoa = 271^3/8. In both cases /l(^) 
is oscillatory, with a wavelength 27t/Qo controlled by C(r), and the lower envelope of the 
minima indicates interface attraction. However, as a consequence of tuning (commensura- 
bility), we find that the attractive minima of ([^ are about one order of magnitude deeper in 
the (111) case. Since the final surface melting behavior is determined by /sr(^) + /l(^), the 
competition between SM and NM reduces, in our picture, to a matter of: a) how large fsR 
(generally repulsive) is; b) how important the liquid surface oscillation 6puv is; and c) how 
bad, or good, is the relative tuning of the two periodicities 271 /Qo and a. The binding energy 
between the SL and LV interfaces is given by A7 = 7sl + 7lv — 7sv = /o°°[/sr(^) + fhi.^)]^^- 
The above result suggests that /l contributes much more to A7 in a well-tuned case. In 
a poorly packed, poorly tuned face, the short-range repulsion cannot easily be reversed. 
In this case, fluctuations are very effective in washing out the shallow layering minima [0, 
whence a negative overall A7, leading to SM. In contrast, the good tuning of a well-packed 
face can cause /l to prevail and yield the final A7 > typical of a NM surface. 

We have tested this qualitative picture by a microscopic ID model calculation. We 
extend the Tarazona- Vicente (TV) scheme — which demonstrates layering of a fluid near 
a wall — to describe a fluid confined between a wall (the LV interface) and a density-wave 
solid with adjustable periodicity. A ID lattice-gas fluid (occupation pi = 0, 1, lattice spacing 
equal to 1), is described by a grand canonical free energy 

n{p,}=F{p,} + J2iV^~p)p^. (2) 
i 

We take 

F{pi} = keT [pi In Pi + (1 - pi) ln(l - p,)] - aptPi+n (3) 

so as to describe a simple fluid. The first term in (^) is entropic, and the second is a 
gradient, representing short-range order in the fluid, with periodicity 2t{/Qo ~ n. The 
external potential Vi consists of a hard wall at i = and of a semi-infinite "crystal" of 
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periodicity a = \n located at some very large distance i = L away. The hard wall generates 
a density oscillation decaying into the liquid from the left with periodicity n. The crystalline 
perturbation generates a second oscillation decaying into the liquid from the right, with 
periodicity Xn, where A is a detuning parameter mimicking the crystalline interplanar spacing 
of a general surface orientation. The values of all the p^'s are determined so as to minimize 
Q. As the crystal is moved closer and closer to the wall by reducing L, the fluid in between 
is modulated with both periodicities, and so is the excess of surface free energy 7(-L) = 
fl{L) + PL, where P = kBTln{l — p) + otp^ is the lattice gas pressure f^. The resulting 
surface free energy difference A7(L) = 7(L) — 7(00) in the case a/ksT = —4 (as used by 
TV for a strongly structured fluid) is shown in fig. ^ together with density profiles. When 
the two periodicities are in tune (a = 2tt/Qo), there is a sequence of attractive minima as 
a function of the thickness L (upper panel). Trapping in the deepest of these minima will 
correspond to NM, as discussed. When tuning is absent (lower panel), the minima are still 
found, but there is no net attraction left. Prevalence of interface repulsion will favor SM in 
such a case. 

The value of these concepts, if correct, should be readily recognizable in experiments 
(not yet available), or at least in a realistic MD simulation. For this purpose, we have 
conducted parallel simulations on Au(lll) and Au(llO), the former being a well known 
NM face [|lOl, and the latter a SM face [jl^]. We have used the "glue model" many-body 
potential [|1^ , which reproduces correctly, among other things, the bulk melting temperature 



(within ~ 1%), the solid- vapor |T2[ and the liquid- vapor |^ surface energies, the main surface 



reconstructions of gold [0, as well as the NM and SM behavior of these two surfaces 
P JT3|JT^ . MD simulations were done on a variety of systems, in the form of A^-layer slabs 
(32 < N < 80) with M atoms per layer (48 < M < 280), one free surface and lateral 
periodic boundary conditions. Special care was taken in determining the effective bulk 
melting temperature, which in a cell-periodic system is both size and orientation dependent 
ig [for example we find T„ = 1355 ± 5 K for a (111) bulk cell with = 80, M = 56, and 
Tm = 1327 ± 3 K for a (110) bulk cell with A^ = 32, M = 280. These values are close to the 



actual Tm of Au, 1336 K]. Near T^, we generate configurations [which can be of equilibrium 
on Au(llO), and near equilibrium on Au(lll)], where a certain number n of surface layers 
are melted. The density profiles of the two interfaces, SL and LV, are shown in fig. |1|. 
The two oscillations, one attached to the LV interface, the other to the crystal, are clearly 
visible. As anticipated, their wavelengths are quite close for Au(lll) {2it/Qo = 2.30±0.05 A, 
a = 2.39 A), and very different for Au(llO) (a = 1.46 A). 

In order to obtain information on the behavior of the total free energy as a function of £, 
we have developed a new method based on non-equilibrium recrystallization runs. Starting 
from a slab configuration with a large number of melted layers, the temperature T is suddenly 
reduced from an initial value Tj > to a final value T < T^. The SL interface moves 
towards its final equilibrium position with a velocity di/dt (fig. ^). As long as the liquid film 
is thick, the interface motion is approximately uniform. At small thickness, however, we find 
a final speedup for the (111) face but, in contrast, a final slowdown for the (110) face. The 
reverse experiment, i.e. fast melting, can also be done, and the reverse behavior is observed: 
initially, the melting front moves out quickly on (110), slowly on (111). By analyzing the 
behavior of £ as a function of time t, we extract the effective SL-LV interaction free energy 
/(£), the average now fully including interface fluctuation effects. The driving force on the 
SL interface is given by —dJ^/di, where dJ^ is the change in the total free energy per unit 
area when the liquid thickness is changed from i to i + di: 

dJ^ = pL{l - T/Tm)d£ + /(£ + dl) - m, (4) 

where p is the density of the liquid, and L the latent heat of melting. Close to Tm, the 
velocity may be assumed to be proportional to the driving force (classical body in a viscous 
medium) : 

df 

- = -A[pL{l-T/T^)+f'{e)] (5) 

For £ > 20A, f'{i) appears to be negligible, and we obtain Am = 3.0 x 10^^A^s~^meV~^ 
and Alio = 4.0 x 10^^A^s~^meV~^. We checked that they do not depend on T in the region 
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of interest (within lOOK from T^). By assuming f{i) = — A7 exp(— ^/^), eq. can be 
easily integrated, giving 

A7 \ 



e{t) =ein ^Ce"*/" 



(6) 



^pL{l-T/TjJ 

where r = C,/[ApL{l — T/Tm)] and C is the integration constant, fixed by the initial liquid 
thickness £(0). In the equilibrium limit t — > +00, and if A7 < 0, eq. (|]) reduces to the 
well-known logarithmic growth law of £ as a function of T 0. From previous equilibrium 
or quasi-equilibrium simulations [^|13],|14| we extract (using also the relationship between 
maximum overheating temperature and A7 [0) A7110 = —8.5 ± 3.0meV/A^, ^no = 2.3 ± 
0.2 A, A7ni = +7.5 ± 2.0meV/A2, ^ni = 11 ± 1 A. Together with p = 0.0562 A'^, L = 
0.112 eV/atom (known values for this potential) we calculate the theoretical i{t) from eq. 
(^). This reproduces simulation data fairly well (fig. ^), except where —dJ-'/di becomes 
small and dynamics slows down [e.g., i ^ 15A for (110)]. Clearly, the overall effective 
interaction free energy is repulsive in the (110) case, and attractive in the (111) case. Its 
magnitude is surprisingly large, as much as 10% of the total surface free energy (71meV/A2 
for liquid Au at Tm [0]). The individual layering oscillations are generally not observable in 
the recrystallization interface velocity, washed out by the large intrinsic interface width, and 
by fluctuations. However, the large value of shows that the (111) attraction is liquid- 
mediated, and must therefore come from layering, as also suggested by fig. ^ Conversely, 
the small value of ^no reflects the destructive interference of the two oscillations. 

In summary, we have argued theoretically, and demonstrated both in a ID model and 
realistic recrystallization simulations the changeable role of the layering on surface melting. 
The strong effect of the liquid surface density oscillations, and their tuning-detuning effect 
has been demonstrated. Through X-ray detection of the magnitude of these oscillations, 
now possible |jl6|, these results should shed new light on the high-temperature behavior of 



solid surfaces. 
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FIGURES 

FIG. 1. Comparison between the (x, y averaged) density profiles of the (111) and (110) surfaces 
of LJ and Au, when wetted by a hquid film of thickness £ ~ 30 A (or ^ ~ 12 LJ units), obtained by 
constant energy MD simulations. The film thickness is stabilized by energy conservation. Layering 
oscillations induced in the liquid film by the solid are present in all cases. The metal also exhibits 
a strong layering on the opposite side of the film, induced by the liquid surface. Its wavelength 
is rather different from that of the solid-induced oscillation in the case of Au(llO), but very close 
for Au(lll). Upon releasing the constant energy constraint at T ~ T^, Au(lll) is unstable and 
recrystallizes [111] . 

FIG. 2. Variation of the surface free energies A7 (with respect to the infinite separation limit) 
as a function of the interface separation L, using the lattice gas models described in the text. 
Upper panel: (lll)-like case (in-tune perturbations). Lower panel: (llO)-like case (out-of-tune 
perturbations). The lower envelope of the two curves shows strong interface attraction for "(111)", 
and weak attraction followed by repulsion at shorter distances for "(110)". Corresponding typical 
plots of densities pt as a function of lattice site are shown in the insets. In the "(110)" case, the 
periodic potential induces crystal-like density oscillations for i > L. 

FIG. 3. Liquid thickness £ as a function of time t as obtained from non-equilibrium MD 
quenching simulations from T = 1500 K to T = 1250 K ~ 0.93 Tm for Au(lll) (filled dots) and 
Au(llO) (open dots). The starting point was in both cases a liquefied slab, more than 100 A thick, 
sitting on top of a few (111) or (110) rigid crystalline layers. Curves corresponding to eq. (|^ 
are shown as a solid [(111)] and dotted [(110)] line. Recrystallization proceeds at constant (but 
slightly orientation-dependent) speed until £ ~ 40 A. For (111), the speed increases in the final 
stage (leading to a crystalline (111) surface), indicating attraction between the SL and the LV 
interface. The opposite occurs on (110), leading to a surface melted state. In this case the system 
reaches the equilibrium thickness near t ~ 1.2 ns. 
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